Spontaneous- Symmetry-Breaking Mechanism of Adiabatic Pumping 
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We consider heterostructures consisting of regions with a continuous symmetry in contact with 
regions wherein the symmetry is spontaneously broken. The low-frequency dynamics of the corre- 
sponding order parameter are shown to induce nonequilibrium transport, a "pumping," out of the 
symmetry-broken regions, which is governed by the generator of the broken-symmetry operator. 
This pumping damps Goldstone-mode excitations and transfers them beyond traditional (static) 
proximity length scales. Our general conclusions are discussed for specific examples of order pa- 
rameters in helimagnets, charge/spin-density waves, superconductors, and ferromagnets. We carry 
out a detailed calculation of such pumping for spiral magnetic orders in helimagnets possessing a 
duality in the representation of its symmetry-broken states. 
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Transport in hybrid correlated-electron systems dis- 
plays a variety of exciting physical phenomena that are 
governed by collective degrees of freedom. Examples 
include supercurrents in Josephson junctions and per- 
sistent spin currents in exchange-coupled ferromagnetic 
spin valves. Interesting effects can also occur when the 
nonequilibrium quasiparticlc transport is initiated and 
manipulated by controlling the collective degrees of free- 
dom. We show in the following that such nonequilibrium 
transport should be generated by any symmetry-breaking 
(SB) order parameter (OP) that is varied by external 
fields. As a specific example of pedagogical value, we 
focus on the spin and momentum pumping induced by 
sliding helimagnets in contact with normal metals. 

Consider an "island" characterized by a local finite or- 
der parameter (f>(r) describing some spontaneous SB em- 
bedded into a "host" with vanishing <p. Electron corre- 
lations maintaining the finite OP in the island can pen- 
etrate into the host giving rise to equilibrium proxim- 
ity effects in hybrid structures. We wish to study the 
nonequilibrium case when the time-dependent OP </>(r, t) 
is driven by an external field. It will be shown that the 
host quasiparticles respond to the order parameter dy- 
namics close to the SB region, generating nonequilibrium 
transport that can propagate beyond the static proxim- 
ity length scales. We develop a general phenomenological 
theory of dynamic coupling between the host and the is- 
land in terms of nonequilibrium pumping generated by 
the OP variation. This pumping generically exerts re- 
action forces that can be important for the OP dynam- 
ics. For demonstration purposes, we consider specific ex- 
amples of helimagnets (HM), charge/spin-density waves 
(CDW/SDW), superconductors (SC), and ferromagnets 
(FM) in contact with normal metals (NM). Our conclu- 
sions may find applications in various areas of physics. 

Spin emission by precessing ferromagnets is already a 
well-studied phenomenon, which has recently attracted 
a renewed interest, leading to many useful insights af- 



ter it was formulated as a parametric-pumping pro- 
cess pj. In Ref. 0, we showed that spin and charge 
pumping by magnetic superconductors reflect their spin- 
pairing symmetry in experimentally-observable quanti- 
ties. Here we show that the mean- field dynamics of any 
symmetry-breaking OP generate nonequilibrium quasi- 
particle transport and propose an approach to simplify 
its computation. We apply the theory to helimagnets and 
charge/spin-density waves, as well as discuss the already 
studied cases of ferromagnets and superconductors in the 
context of a generalized formalism. To this end, let us 
consider an arbitrary SB island with partition function 
expressed as a functional integral over the OP field 
Z = J T>[<p(r)] exp[-/3{.F[0(r)] + / dr<p(r)H(r)}], where 
(3 is the inverse temperature, J- is an effective free-energy 
functional, and H(r) is an external field coupled to the 
OP. Suppose the free energy T is invariant under a con- 
tinuous transformation group U(6) parametrized by 9 
[U(0) = 1], J-[U<j>] = J- [(/)], defining a symmetry that 
is broken by the states with nonvanishing OP <p. In the 
absence of the external field H, U applied to a reference 
ground state (f>o generates a manifold $ of states min- 
imizing T[cj)]: 4>0 (r) = U(0)<j)o(r). We are interested in 
the symmetry subgroup with nontrivial [apart from U(0)] 
and unique transformations on 4>q: U{6)4>q — U(6')<po 
only and only if 6 = 6'. A slowly- varying H(r,t) can 
adiabatically move </>[W(r)] through giving us a han- 
dle to control the dynamics of the OP. Alternatively, a 
small oscillating component of Ti on top of a larger static 
component 7io can set off a resonant motion of <f> near 
<f>[Ho(r)], as in the case of the ferromagnetic resonance. 

The essential ingredient in our model is an effective 
mean- field Hamiltonian experienced by the host, which 
is modulated by the time-dependent OP cf>(t) in the is- 
land. We view the SB island as a scatterer perturbing 
the state of the host. Coupling to the host also affects 
its properties, which should be taken into account self 
consistently. We wish to consider the host response to a 
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time-dependent Hamiltonian H[(j>(r,t)]. To be specific, 
we consider an effective Hamiltonian H = Ho + H' de- 
scribing electrons in a normal Fermi liquid, which are 
perturbed by the presence of an SB island. H is an un- 
perturbed Hamiltonian and H'[4>(r)] is a coupling to the 
SB island. At sufficiently low temperatures, the island is 
characterized by a given position <f)g (r) on the manifold $ 
of the ground states, which is controlled by the small ex- 
ternal field h and is labeled by 9 defining a transformation 
from a certain reference state </>o(r). The Hamiltonian is 
also labeled by 9 correspondingly: H'(6) = H'[(j>g(r)], 
with each 9 assumed to define a unique H' . Suppose 
H'(9) can be generated from H'(0) by a certain unitary 
^-dependent transformation on the host degrees of free- 
dom, and Hq is invariant under this transformation, 

H'{9) = U(9)H'(0)U\9) , [H , U(0)} = 0, (f ) 

where we naturally use the same symbol U . 

To give more substance to our discussion, we consider 
several concrete examples shown in Fig. The Hamil- 
tonian for host electrons scattering off the SB island and 
locally experiencing its OP is given in these examples by 
H ' = J dr over the integrands 

(a) ^(r^ryQ-r, (b) *t $ a (r)C(r)e«' r , 

(c) *{(r)*t(r)A(r) , (d) ¥ a * a ,(r) [* aa , • e xc (r)] (2) 

plus Hermitian conjugates in (a), (b), and (c), and the 
sum is implied over repeated indices. Here ^c's are 
spin-a electron field operators, h(r), C(r), and A(r) are 
complex-valued HM, CDW, and SC gap functions, re- 
spectively, e xc is the local ferromagnetic exchange field, 
and er is a vector of the Pauli matrices. We take C(r) to 
be uniform along the direction of Q. The SDW Hamil- 
tonian can be obtained from that of the CDW by adding 
an extra factor (— 1) Q under the sum. The CDW can 
result from a Peierls instability in quasi-one-dimensional 
electron gases, leading to Q = Ikp in terms of the Fermi 
wave vector fcj? Q. The OP's, which self-consistently 
define h, C, A, and e xc experienced by quasiparticlcs, 
are given respectively by the HM spiral order ('I'j'I'j) 
0, the CDW condensation amplitude (^t^a) [!■ the 
Cooper-pair amplitude (^I , |\E'|), and the itinerant mag- 
netization (^(T att <$ a '}. We assume that the ferromag- 
netic exchange field can point in arbitrary direction m, 
e xc (r, t) = e xc (r)m(t), breaking the spin-rotational in- 
variance of the Hamiltonian H, and the gap functions 
can have arbitrary phase ip, /(r, t) = f(r)e~ lv ^\ break- 
ing the spin-rotational invariance around the z axis for 
the helimagnct (/ = h), translational invariance for the 
charge/spin-density wave along Q (/ = C), and gauge 
invariance for the superconductor (/ = A). Helimagnets 
whose gap function h(r) is uniform along Q will be of a 
special interest for our discussion because of the dual rep- 
resentation of their symmetry-broken states: Any state 



determined by an overall phase of h(r) can be obtained 
from a reference state with zero phase by either a rotation 
in spin space around the z axis or a spatial translation 
along the direction of Q, see the upper inset of Fig. [3 
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FIG. 1: Schematics of SB heterostructures. (a) HM dynamics 
can be viewed either as a rotating or sliding spiral order with 
frequency lj = dt<p and wavevector 2tt/Q. These dynamics 
induce both momentum and spin flow into the normal metal, 
(b) A sliding CDW felt by an NM host that locks to its motion. 
CDW phase variation translates into velocity v = dt<p/Q. 
The dynamics thus pump momentum, (c) Charge current 
I c pumped by the varying pairing-gap phase tp(t) induces a 
charge imbalance in the normal metal corresponding to the 
voltage V = (h/2e)dt<p. (d) Spin pumping I s by the rotating 
magnetization direction m(t) establishes spin splitting e = 
h\m x dtm\ along the rotation axis. 

The mean-field OP dynamics translate into the time 
dependence of H'. Let us write the corresponding trans- 
formation along a given path parametrized by a single 
parameter 9(t) — uit, where u> is the velocity along this 
path, as U(t) = exp(-iGuit), in terms of the genera- 
tor G = ideU(9)\e = Q. Due to adiabaticity, any general 
motion can instantaneously be described by a constant- 
speed trajectory. The time dependence of H = Ho + H' 
is formally eliminated by applying the inverse transfor- 
mation, U*{t), to the host degrees of freedom and corre- 
spondingly adding a new term to the Hamiltonian, de- 
termined by the broken-symmetry generator: 

H{t) -> H(0) - hwG, (3) 

subject to Eqs. JIJ. This can be loosely interpreted as 
follows: The effective Hamiltonian experienced by the 
host medium due to the dynamic SB island is time inde- 
pendent in the "reference frame" moving with the island. 
The additional term hwG in Eq. reflects the inertial 
forces in the moving frame of reference and is given in 
our toy models, Eqs. respectively by 

(a) S z dw or (P • Q/Q 2 )c^, (b) (P • Q/Q 2 )d^, 
(c) (hN/2)d tV > , (d) S • (m x d t m) , (4) 



3 



where S is the total-spin operator, P is the total- 
momentum operator, and N is the total-particle- number 
operator. The example of spiral order is special in that 
the time dependence of the Hamiltonian J2t) can be re- 
moved by two different transformations: either rotation 
in spin space or translation in real space [if Q- d T h(r)=0]. 
The translation generator |0p) can be removed by a 
gauge transformation (assuming the usual nonrelativis- 
tic Hamiltonian with kinetic energy p 2 /2m, per particle 
with momentum p, and electron-electron interactions), 
which is consistent with the Galilean invariance dictat- 
ing that there should be no inertial forces in a reference 
frame moving linearly at a constant speed. 

If U is a global-symmetry operator for the entire host 
medium, the new term in the Hamiltonian © leads to a 
stationary state that could be different from that of Hq 
even very far from the SB island. The simplest example 
is the CDW: The time dependence of the gap-function 
phase ip(t) = Lot corresponds to sliding of the effective po- 
tential experienced by the electrons at velocity v = lo/Q 
along the direction given by Q, see Fig.^b). Due to the 
gauge invariance, eigenstates of the transformed Hamil- 
tonian H(0) — vPq, see Eq. (@p), can be obtained from 
the eigenstates of H(Q) by shifting the momenta of all 
electrons by k = mv/H along Q. Transforming them 
back into the "laboratory frame of reference" shifts the 
wavefunctions in space without affecting the momenta, 
and we thus see that the entire host is sliding with veloc- 
ity v locked with the CDW, as can be naively expected. 
If the CDW starts slowly accelerating in the rest frame, 
the host will acquire momentum by "pumping" of finite- 
momentum electron-hole pairs from the CDW conden- 
sate that originates near the CDW region and propa- 
gates across the sample. An SDW can similarly pump 
momentum-carrying electron-hole pairs, an SC pumps 
Cooper pairs, an FM pumps spin-polarized electron- 
hole pairs, and an HM pumps both spin-polarized and 
momentum-carrying electron-hole pairs, as explained in 
the following. 

Beforeproceeding, we point out that the spin-pumping 
physics |l| can be recast as merely another example of the 
SB mean-field pumping: According to Eq. (@p) , magnetic 
direction m(t) rotating with frequency to spin polarizes 
the nonmagnetic structure in contact with the ferromag- 
net. In steady state, the nonmagnetic region is spin split 
along the rotation axis by energy fko. This stationary 
spin accumulation, in turn, interacts with the FM is- 
land inducing spin currents from the NM region, which 
have to cancel the pumped spin currents Q. Knowing 
a steady-state solution and calculating the dc transport 
corresponding to the nonequilibrium state of the system 
thus allows one to find the pumping circumventing the 
time-dependent transport problem. Another simple ex- 
ample that was already considered in literature [2( is the 
superconducting island in Ohmic contact with a normal 
metal, see Fig. ^d). The time-dependent gap function 



A(0)e _4V '*- ) effectively induces a voltage imbalance across 
the junction, determined by the Josephson equation. 

It is instructive to discuss in some detail the HM|NM 
dynamic system. The upper inset of Fig.|3depicts the he- 
limagnet with h(r,t) — h 5(z)e~ luJt and Q = Qx. Since 
the Hamiltonian J2J}) is translationally invariant along 
the y direction, we suppress this axis and consider a two- 
dimensional transport problem in the xz plane (with the 
HM magnetization lying in the xy plane). When such a 
helimagnet is treated as a perturbation embedded into a 
free-electron gas, the spiral dynamics have the following 
stationary solutions corresponding to Eqs. J?p,) : a Fermi 
sea spin polarized along the z axis by energy splitting of 
fko or a Fermi sea drifting with velocity lo/Q along the 
x axis, see lower inset of Fig. (or a linear superposi- 
tion of the two: e.g., spin-polarized by ftio/2 and drifting 
with velocity lo/2Q). Using Fermi's Golden Rule, it is 
straightforward to calculate the backflow of the corre- 
sponding spin accumulation or momentum imbalance to 
obtain the nonequilibrium pumping. We compute that 
the spin-imbalance scattering must be compensated by 
the spin and momentum pumping (at zero temperature, 
per HM unit length) 

dS z /dt = LoS p V{q) and dP x /dt = LoS p V{q)Q , (5) 

where V{q) = J™ 8 ^ 9 " 1 ) d{/y/l - {cos( - q) 2 for < 
q < 2 and V(q) = for q > 2, in terms of q — \Q/kp\, and 
the integration is performed over the angle of incidence 
£ with respect to the x axis. S p = \ho\ 2 m 2 /(TT 2 H 2 kp), 
where m is the electron mass and kp is the Fermi 
wavevector. [It is straightforward to generalize Eqs. © 
to a layered geometry in three dimensions by integrat- 
ing over the phase space along the y axis.] We plot the 
pumping strength V dependence on Q (which can be ex- 
pressed in terms of elliptic integrals) in the main panel 
of Fig. In the long-wavelength limit, Q — > 0, when 
the HM becomes a uniform (5-function FM, V diverges 
logarithmically. Mathematically, this is anticipated since 
/ d£/sin£ = ln[tan(£/2)] diverges at C —> 0,7r. Physi- 
cally, it means that spin pumping by a 5-function ferro- 
magnet is not a smooth function of the exchange field 
h at h — 0, so that the perturbation expansion breaks 
down. This problem has been already encountered in lit- 
erature and can be dealt with by properly including 
a cutoff for electrons that are incident on the magnetic 
layer at a shallow angle ( ~ hokF /2ep. For short wave- 
lengths, Q — > 2kp (but Q < 2kp), pumping approaches 
a finite value V = tt/2, although the integration range 
shrinks to zero. At Q > 2kp, the Golden- Rule integra- 
tion range vanishes and so does pumping, exhibiting a 
step at Q = 2fcp. For a three-dimensional system, how- 
ever, pumping vanishes at Q — > 2kp and there is no step 
in V, since only the largest "slice" of the Fermi gas in 
the k x k z plane with k y — gives a finite contribution. 

For pumping corresponding to the drifting station- 
ary solution, we obtain Eqs. |J5J but with the integral 
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Q/2k F 

FIG. 2: Main panel: V defining the HM pumping, Eq. Q, 
as a function of Q/2kF (with the y axis suppressed in the 
calculation). Upper inset: Spiral magnetic order in the heli- 
magnet. Arrows show the local exchange field in the xy plane 
for h(r) — ho8(z)e~ ltp and Q = Qx, where ho is real-valued 
and Q > 0. Phase tp fixes the HM position in the x direc- 
tion, ip can be varied by either spatial translation along the 
x axis or spin rotation around the z axis normal to the xy 
plane. Lower inset: Two distinct stationary states of the nor- 
mal metal corresponding to the HM dynamics. Left is the 
spin-polarized Fermi gas with energy splitting e = huj along 
z and right is the Fermi gas drifting along x with velocity 
v = u)/Q. 

P(q) = {21 q) / O co '" 1(, - 1) dC cos C/ VI - (cos C - q) 2 for 
< q < 2, which at first appears to be different from 
that in Eqs. JSJ. It turns out, however, that the integrals 
are numerically indistinguishable (although we were not 
able to prove it analytically). The pumped flows thus 
do not depend on the choice of a stationary solution, as 
expected to lowest order in frequency. It is however not 
clear which of the various stationary states is established 
in practice. Unlike pumping, this depends on the proper- 
ties of the normal metal away from the HM|NM contact. 
If, for example, the electron gas cannot drift with the heli- 
magnet because of boundaries or disorder scattering, HM 
dynamics will induce the spin accumulation Huj along the 
z axis, which must cancel the nonequilibrium spin pump- 
ing, in the absence of spin-relaxation processes. If, on the 
other hand, NM electron spins relax without momentum 
scattering, the electron gas will drift with velocity oj/Q 
along the x axis. 

HM dynamics assumed in our theory can be excited 
by a magnetic field h z along the z axis. In the case of an 
isolated dissipationless helimagnet, h z induces HM dy- 
namics dfip = jh z (t), where 7 is minus the gyromagnetic 
ratio. In order to see that these dynamics dissipate en- 
ergy in the presence of an NM reservoir, we calculate the 
imaginary part of the magnetic response function along 
the z axis: 

Im X H = 7 2 ^/ [tfS p K) 2 + oj 2 ] , (6) 



where K is the easy-plane anisotropy constant, assuming 
the HM prefers to lie in the xy plane without an easy axis. 
Energy dissipation is an experimentally-observable man- 
ifestation of the pumping. Consequently, this demon- 
strates that the framework presented in this article can 
lead to easy calculation of experimentally-relevant quan- 
tities. Note that the spin pumping by a spiraling HM 
is different from the persistent spin currents predicted 
in Ref. 0] for an equilibrated HM in the presence of a 
uniform magnetic field lying in the magnetization plane. 

Let us now summarize a general approach for cal- 
culating the adiabatic pumping by any time-dependent 
Hamiltonian breaking a symmetry of the system: 1) The 
broken-symmetry operator U satisfying Eqs. Q is iden- 
tified, 2) the system is transformed into the moving ref- 
erence frame and solved for the stationary state after 
adding the new term in the Hamiltonian, which is pro- 
portional to the generator of U, Eq. (J3J, 3) the station- 
ary nonequilibrium configuration induced by the addi- 
tional term in the transformed Hamiltonian, and the cor- 
responding backflow transport, are calculated in the lab- 
oratory frame. As this transport must be canceled by the 
pumping in the steady state, the procedure automatically 
gives (minus) the pumping due to the time-dependent 
Hamiltonian H'{t). This assumes that the total trans- 
port is given by the sum of the pumped current gener- 
ated by the time-dependent OP and the backflow current 
induced by the nonequilibrium build up, which requires 
adiabaticity of the dynamics. We remark that it is not 
necessary that the steady state is reached in practice, 
since the pumping close to the SB island is not affected 
by the processes deep in the host. For a periodic varia- 
tion of H'(t), the average pumping depends on the tra- 
jectory in the spontaneous-SB manifold. In particular, if 
the trajectory can be parametrized by a single periodic 
real-valued parameter, the net pumping over a cycle van- 
ishes 0- In all of our examples, a constant-rate phase 
winding, tp(t) = tot, (or a steady magnetization rotation) 
can be described by two out-of-phase periodic parame- 
ters, cos(ip) and sin(ip), producing a finite net pumping. 

Finally, we wish to comment on the reaction of the 
pumping on the OP dynamics. Since the process dis- 
sipates energy, the reaction of the medium on the SB 
island should damp the Goldstone modes restoring the 
broken symmetry. This damping can be calculated by 
first finding the pumped transport and then imposing the 
conservation laws in the classical equations of motion of 
the OP, as we did above for the helimagnets. A related 
phenomenon is the dynamic interaction between differ- 
ent SB regions: The nonequilibrium pumping induced by 
the OP dynamics can propagate from one SB island to 
another, damping Goldstone modes in one and exciting 
them in the other. Such dynamic interaction can be long 
ranged in comparison with the static interaction that re- 
quires quantum proximity. This formalism have already 
led to a considerable advancement in understanding of 
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the ferromagnetic dynamics in hybrid structures [lj . 
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